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Abstract.  In  the  present  paper  we  offer  a model  of  the  growth  kinetics  of  Sii_.vGe,  layers  from 
silane  and  germane  molecular  flows.  The  model  considers  disintegration  of  silane  and  germane 
molecules  from  chemisorption  to  formation  of  Si  and  Ge  adatoms  and  their  embedding  in  a 
crystal  lattice.  The  numerical  modeling  of  the  epitaxial  process  has  shown  good  agreement  with 
the  experimental  data  and  has  allowed  to  explain  peculiarities  of  the  experimental  dependencies. 
The  effective  frequencies  of  SiH3  and  SiH2  disintegration  for  the  considered  epitaxial  method 
were  estimated. 


Introduction 

In  the  last  decade  the  chemical  vapour  deposition  from  silicon  and  germanium  hydrides 
at  low  pressures  has  been  actively  used  for  Sii_.vGe,  layer  epitaxy.  Despite  the  large 
number  of  experimental  works  the  number  of  theoretical  researches  is  insignificant,  and 
the  available  models  describe  experimental  data  inadequately.  Therefore,  as  remarked 
by  L.  T.  Vinh  with  co-workers  “improved  models  thus  seem  to  be  indispensable”  [1], 
The  not  good  enough  agreement  of  the  theory  and  experiment  is  possible  connected 
to  an  extremely  simplified  consideration  of  chemical  processes  on  a growing  surface. 
On  the  other  hand,  a detailed  description  of  the  kinetics  of  surface  reactions  inevitably 
will  result  in  an  increase  of  a number  of  parameters  describing  the  rates  of  elementary 
processes.  The  task  is  complicated  also  by  that  the  values  of  the  majority  of  rate 
constants  disintegration  of  SiH,,  and  GeH„  molecules,  unfortunately,  are  unknown. 

In  the  present  paper  we  offer  a model  of  the  growth  kinetics  of  Sij  _xGex  layers  from 
silane  and  germane  molecular  flows.  The  model  considers  disintegration  of  silane  and 
germane  molecules  from  chemisorption  to  formation  of  Si  and  Ge  adatoms  and  their 
embedding  in  a crystal  lattice.  Using  the  published  experimental  data,  on  the  basis  of 
the  model  we  have  made  estimations  of  rate  constants  of  disintegration  of  SiH3  and 
SiH2  molecules. 

1 Model 

A comprehensive  analysis  of  the  process  of  disintegration  of  SiH„  molecules  on  a surface 
Si(100)  is  made  by  S.  M.  Gates  with  co-workers  [2].  They  found  out  that  the  complete 
mechanism  of  silane  disintegration  is  reactions: 


sSiH4  MH,  MH2 

SiH4  + 2 -»  SiH3  + H (1),  SiH3  + 1 ->  SiH2  + H (2),  2SiH2  — > 2SiH  + H2  (3), 

„P> 

MH,  , , "h  , , , 

SiH2  — » SiH  + H 4),  SiH  4 Si  + ±H2  (5). 
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Taking  into  account  a chemical  similarity  of  SiH4  and  GeH4  molecules  it  is  possible 
to  expect  that  the  disintegration  of  germane  will  occur  by  an  analogous  way  [3],  If  on 
a growing  surface  there  are  SiH2  and  GeH2  molecules,  these  molecules  possibly  can 
react  by  the  reaction: 


SiH2  + GeH2 


SH  + GeH  + H2  (6). 


An  appropriate  system  of  kinetic  equations  for  dimensionless  average  concentrations 
of  adsorbed  particles  is  (M  = Si,  Ge): 


ae, 


dt  ~ 1 n, Unbi  ~ %H5yMHTv> 

~di~L  = ~ ~~  2 ,yMH , %iH 3 9G c It : , 

de 


or1  = 4Vll/2MH,  + + 2^MH,0SiH20GeH2  - 


de 

~ai 


= E (2  M\MH4^  + 2i/MH3 %H3 0>'  + 2i/MH, 00  % (7) 

M=Si,Ge x 7 


- rivt^M  “ rdM> 

0V  + 0Mh3  + #MH,  + %H  + 

Ov  — @bl  ”t“  @nbU 
dbi  = 3 (^SiHi  + 0GeH3)  • 


The  system  (7)  is  based  on  the  following  considerations,  (a)  Monomolecular  ad- 
sorption takes  place  on  a growing  layer  surface,  (b)  Lateral  interactions  between  atoms 
and  molecules  on  a growth  surface  are  absent,  (c)  Chemisorption  of  silane  and  ger- 
mane molecules  requires  two  vacant  not  blocked  adsorption  sites,  (d)  The  GeH3  and 
SiH3  molecules  block  3 vacant  adsorption  sites,  (e)  The  desorption  rate  of  hydrogen 
atoms  from  an  alloy  surface  linearly  depends  on  Ge  content,  (f)  The  frequency  of  atom 
embedding  into  a lattice  (r)  depends  on  concentration  of  atoms  Si  and  Ge  as  described 
in  [4],  Growth  rate  and  Ge  content  of  an  alloy  layer  can  be  calculated  by  the  formulas 


9Gc 

®Ge  + %i  ’ 


2 Results  and  Discussion 


Because  the  silane  properties  are  investigated  in  more  detail,  we  have  focused  the  efforts 
on  modeling  of  a growth  from  SiH4.  In  order  to  define  values  of  kinetic  constants  the 
experimental  dependences  of  growth  rate  of  layers  on  temperature  [1,  5]  and  TPD 
spectra  of  hydrogen  [2]  have  been  used.  As  is  visible  from  Fig.  1,  the  dependence  V(T) 
has  a kink  at  the  temperature  about  620-630  °C.  The  V{T)  slope  above  this  temperature 
is  small  and  the  activation  energy  approximately  equals  0.1  eV.  Such  a value  of  the 
activation  energy  correlates  well  with  the  activation  energy  of  disintegration  of  SiH3 
molecules  [2],  Thus,  the  increase  of  the  growth  rate  at  T > 620— 630  °C  is  limited  by 
the  disintegration  of  SiH3  molecules.  The  fitting  of  the  theoretical  dependencies  V(T) 
to  experimental  data  has  allowed  to  determine  the  disintegration  frequency  of  SiH3 
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Fig  1.  Experimental  (curve)  and  theoretical  (triangles  — [1],  circles  — [5])  growth-temperature 
dependencies  of  the  growth  rate. 


Surface  temperature  (K) 

Fig  2.  TPD  data  showing  hydrogen  desorption  for  silane  adsorbed  on  the  Si(100)  [2], 


molecules  r-'MH-,  — 2.5  x 10°e_0087(eV^ir(s_1).  This  result  coincides  with  conclusions 
of  S.  M.  Gates  with  co-workers  [2], 

The  definition  of  disintegration  frequencies  of  SiH2  molecules  is  a more  difficult  task. 
The  difficulty  consists  in  the  fact  that  the  disintegration  of  these  molecules  proceeds 
by  two  ways  simultaneously.  In  order  to  determine  these  parameters,  we  have  used 
TPD  spectra  of  hydrogen  [2],  As  is  visible  from  Figs.  2 and  3,  these  spectra  have  two 
maxima  at  700  and  800  °C.  The  occurrence  of  the  first  desorption  maximum  is  caused 
by  reaction  between  two  SiH2  molecules,  the  second  maximum  is  connected  with  the 
simple  hydrogen  desorption  from  the  Si(100)  surface.  Fitting  the  first  peak  we  have 
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Fig  3.  Calculated  hydrogen  TPD  spectrum  for  silane  adsorbed  on  the  Si(100). 

found  ~ 5 x \QX1  e~2  5('eV')/kT (s~x).  The  fitting  of  the  second  peak  has  allowed 
to  specify  the  frequency  of  hydrogen  desorption  from  the  Si  (100)  surface  (;/j|  ~ 
8 x 1012e-213(eV^/cr(s-1)).  The  second  peak  has  one  feature:  the  low  temperature  half 
of  the  peak  is  narrower  than  the  high-temperature  one.  The  extension  of  this  peak  in  the 
high  temperature  field  is  apparently  connected  with  the  monomolecular  reaction  (4). 
For  the  kinetic  constant  of  this  reaction  we  have  found  j/j ^ ~ 5 x 102e-0'6(eV}/*r(s-1 ). 
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